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ABSTRACT 

For high-precision cosmic microwave background (CMB) experiments, contamination from extragalac- 
tic point sources is a major concern. It is therefore useful to be able to detect and discard point source 
contaminated pixels using the map itself. We show that the sensitivity with which this can be done 
can often be greatly improved (by factors between 2.5 and 18 for the upcoming Planck mission) by a 
customized hi-pass filtering that suppresses fluctuations due to CMB and diffuse galactic foregrounds. 
This means that point source contamination will not severely degrade the cleanest Planck channels unless 
current source count estimates are off by an order of magnitude. A catalog of around 40,000 far infra-red 
sources at 857 GHz may be a useful by-product of Planck. 
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1. INTRODUCTION 

If future CMB experiments are to produce high- 
precision measurements of cosmological parameters (Jung- 
man et al. 1996; Bond et al. 1997; Zaldarriaga et al. 1997), 
they must remove foreground contamination from galactic 
dust, synchrotron and free- free emission as well as extra- 
galactic point sources with comparable accuracy (see e.g. 
Brandt et al. 1994; Tegmark & Efstathiou 1996 - hereafter 
TE96; Bersanelli et al. 1996). Fortunately, foregrounds 
differ from CMB fluctuations in several ways, all of which 
can be used as weapons against them, in combination: 

i. Their non-Gaussian behavior can be used to discard 
severely contaminated regions (e.g., bright point 
sources, the Galactic plane). 

ii. Their frequency dependence can be used to subtract 
them out by taking linear combinations of maps at 
different frequencies. 

iii. Their power spectra can be compared with the one 
measured, to fit out remaining foregrounds. 

Extragalactic point sources are one of the most menac- 
ing foregrounds, for several reasons: 

• Their spectral index varies much more than for other 
foregrounds. 

• Many sources exhibit substantial time variability. 

• Their abundance is very poorly known over much of 
the CMB frequency range. 

The variation of their frequency dependence in a random 
way across the sky (from source to source) substantially 
degrades the effectiveness of weapon (ii) (Tegmark 1998), 
and time- variability of both flux and spectral shape (e.g., 
Gutierrez de la Cruz et al. 1998) further complicates sub- 
traction attempts. As a complement to (ii), it is therefore 
important to make as much use as possible of 1. This is 
the purpose of the present Letter. 

If point sources give the only non-Gaussian contribution 
to a map, then a useful way to implement 1 is to discard 



all pixels whose temperature lies more than v standard 
deviations a above the mean. Chosing v — 5 ensures that 
this will falsely reject only a negligible fraction of order 
3 x 10 -7 of all uncontaminated pixels. The accuracy with 
which this method will be able to clean the CMB maps 
from the upcoming Planck satellite (Bersanelli et al 1996) 
has recently been estimated by Gawiser & Smoot (1997), 
Toffolatti et al (1998, hereafter To98) and Guiderdoni et 
al. (1998). Similar estimates have been made for the MAP 
experiment (Refregier et al. 1998). For such high signal- 
to-noise experiments, the rms pixel fluctuations a are not 
dominated by detector noise but by CMB (or, in some 
channels, by galactic foregrounds). We will show that 
much of these fluctuations (which hamper our ability to 
detect point sources) can be removed by an appropriately 
chosen band-pass filter, and that consequently, the thresh- 
old va at which point sources can be removed can be sub- 
stantially lowered. We derive our method in §2, apply it 
to Planck in §3 and discuss our results in §4. 

2. METHOD 

If there are point sources with fluxes Si at sky positions 
given by unit vectors r, , then the resulting sky tempera- 
ture x is 



x(r)=cJ2SiS(i 



,r) + ^a for 



Yi m (r), 



(1) 



where 5 is a Dirac delta function and the spherical har- 
monic coefficients ai rn contain the combined contributions 
from CMB, galactic foregrounds and detector noise. Here 
c is the conversion factor between surface brightness and 
temperature, given by (see, e.g., equation (3) in TE96) 
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where rj = hv jkT cm b w v/56.8 GHz. The observed map is 
x(r) convolved with the beam function B. To maximize 
our sensitivity to the point sources, we wish to convolve 
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it with an additional function W, giving a filtered map 
y(r) = (W -k B * x)(r) (★ denotes convolution). Equa- 
tion (|l|) gives 

y(r) = c Si (W * B) (r t ■ r) + ^ W e B e a im Y em (r ) . (3) 

i Irtx 

Here we have taken both the filter W and the beam B to 
be spherically symmetric, so that they are given by the 
coefficients Wg and Bg in a Legendre polynomial expan- 
sion: W(cos0) = J2e ( 2 I7 i ) W e P t (cos9) and B(cos6) = 
T,e i 2 ^-) B e P e (cos9). Equation (§) tells us that if the 
contribution from overlap of nearby sources is negligible 
(we will see that this is a good approximation for Planck at 
the attainable flux threshold), then the point source con- 
tribution in the direction of the i th source is y(?i) = ASt, 
where the normalization constant A = (W * B)(l)c, since 
T* • r< = 1. This constant can be re- written as 



A = (W*B)(l)c = cY^ 
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In other words, equation (Q) tells us that the peak bright- 
ness of a source in the normalized map y(r)/A gives its 
strength S in flux units (Janskys). This is true for any 
choice of our filter W, so we simply wish to choose W so 
that it minimizes the variance a 2 in this map. Equiv- 
alently, we want to maximize the signal-to-noise ratio, 
which is the ratio of the peak signal of a source after fil- 
tering to the rms fluctuation level a in this region not 
due to the source. These fluctuations come from the ag m 
in equation (0), i.e., from CMB, pixel noise and galac- 
tic foregrounds, all which act as unwanted noise now 
that point source detection is our objective. Modeling 
these as Gaussian random fields as in TE96, we have 
{a* lm at' m ') = Su'S m rn'Cj ot , where C| ot is the sum of the 
power spectra of the CMB, the noise and various Galac- 
tic foreground components. Using this, equation (||) and 
equation (0) , we find that the variance of our point source 
map y(r)/A is 
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It is easy to show that this variance is minimized if we 
choose the filter to be 
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In summary, the outlier removal method can eliminate all 
point sources with flux S > S c = va, where a is given 
by equation ((?]), and v can be chosen depending on the 
desired confidence level of source detection. 

We note that our assumption of symmetric beams is by 
no means necessary. Indeed, our filtering method is read- 
ily generalized to the case of arbitrary beam shape and 
arbitrary known profiles of (resolved) sources, i.e., galaxy 
clusters, in a manner analogous to Appendix A of Haehnclt 
& Tegmark (1996). 



3. APPLICATION TO PLANCK 

We will now illustrate our method with an application 
to the upcoming Planck mission, to see by what factor it 
reduces a. 



3.1. Assumptions about foregrounds, noise and CMB 

For the Galactic foreground power spectra, we use the 
model C t {u) oc c 2 B{iy) 2 £~ 3 (TE96; Bersanelli et al. 1996) 
for dust, synchrotron and free- free emission, where B is 
the sky brightness measured in MJy/sr. We model the 
frequency dependence as B(y) oc jy -015 for free-free emis- 



sion and as B(y) 



oc v 
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for synchrotron emission (see 



Platania et al. 1998 and references therein). For dust, we 



assume B(v) 



3+/3 U e hv/kT 



1), with a dust temper- 



ature T = 20K and an emissivity /3 = 1.5 (Kogut et al. 
1996, hereafter K96). 

The combined DIRBE and IRAS dust maps suggest a 
slightly shallower slope £~ 2 - 5 (Schlegel et al. 1998), but a 
recent analysis of the DIRBE maps has shown no evidence 
of a departure from an £~ 3 power law (Wright 1998) for 
£ < 300, and we will see that only the behavior at low £ 
matters for the present analysis. 

We normalize the power spectra based on the DIRBE- 
DMR cross-correlation analysis of K96, which gives rms 
fluctuations of 2.7/iK for dust and 7.1/iK for free-free 
emission at 53 GHz on the COBE angular scale, in good 
agreement with the DIRBE cross-correlation results for 
the Saskatoon map (de Oliveira-Costa et al. 1997) and 19 
GHz map (de Oliveira-Costa et al. 1998). We normalize 
the synchrotron model to give 11/jK at 31 GHz, the K96 
upper limit, since cross-correlation between the 408 MHz 
and 19 GHz emission indicates a value in this range (de 
Oliveira-Costa et al. 1998). For radio and infrared point 
sources, we use the source count model of To98. 

The noise power spectrum is C% mse = {FWBMa n ) 2 /B 2 
(Knox et al. 1995; TE96), where a n is the rms noise 
in a pixel of area FWHM 2 . We assume Gaussian 
beams of rms width 6, which corresponds to B( = 
e e 2 i(£+i)/2^ where e = FWHM/ (8 In 2) !/ 2 . The values of 
FWHM and a n for the Planck channels are taken from 



http:// astro, estec. esa.nl/Planck/ . 

We compute the CMB power spectrum with CMBFAST 
(Seljak & Zaldarriaga 1996) for three COBE-normalized 
models. We use a standard cold dark matter model 
(SCDM) with n = 1, h = 0.5, h 2 fl b = 0.015 as well as two 
models that agree better with observational data (Wang 
et al. 1998). These are ACDM, a flat model with a cos- 
mological constant (as SCDM except that fl\ = 0.5), and 
OCDM, an open model (as SCDM except that £1 A = 0, 
Q = 0.5, h = 0.65). 



3.2. Results 

The optimal filter W for Planck channel 5 in the ACDM 
model is shown in real space in Figure 1 and in the Fourier 
(multipole) domain in Figure 2. The ACDM results for all 
channels are summarized in Table 1, and it is seen that 
the sensitivity a is improved by a factor ranging from 2.5 
to 18. This can be qualitatively understood from Figure 2, 
which shows that W is essentially a band-pass filter, tar- 
geting those multipoles where the combined fluctuations 
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from CMB, foregrounds and noise are minimal. The table 
shows that the dominant sky signal (CMB for channels 
1-7, dust for channels 8-10) is suppressed by even larger 
factors, at the cost of higher detector noise. The filter- 
ing technique therefore helps the most for high resolution 
experiments where the signal-to-noise level per resolution 
element is substantially greater than unity. This should be 
attainable for instance for several of the upcoming inter- 
ferometer experiments; see Smoot (1997) for a recent re- 
view. Conversely, the filtering helps only marginally (with 
gains below a factor of 2) for the upcoming MAP satellite, 
for w hich we repeated our anal ysis using the specifications 



from http://map.gsfc. nasa.goi\ This is because the MAP 



detector noise is not too far below the fluctuation levels 
from CMB and foregrounds to start with, leaving less room 
than in the Planck case for filtering to improve the situa- 
tion. The same conclusion is drawn in the MAP analysis 
of Refregier et al. (1998). 

Table 1 shows that the number of sources removed 
changes by an even greater factor than a does. This is be- 
cause the differential source counts of To98 are quite steep. 
Channel 10 is especially notable: here Planck should be 
able to produce a catalog with around 40,000 sources, un- 
doubtedly useful in its own right, as compared to a measly 
300 without filtering. Our filtering technique should also 
be useful for constructing a catalog of SZ-clusters from the 
217 GHz map, as described by Aghanim et al. (1997). 

The SCDM and OCDM models give results quite similar 
to Table 1, with gain factors differing by less than 20%. 




FIG. 1 — The convolution filter W is plotted in real space for 
channel 5 of Planck. 

4. DISCUSSION 

We have presented a method for point source removal 
that consists of the following steps: Convolve the CMB 
map with a band-pass filter such as the one in Figure 1. 
Compute the resulting rms fluctuation level a and inter- 
pret all positive fluctuations exceeding va as point sources. 
Revert to the original map and discard the contaminated 
pixels (within a beam size or two of each source) . 

The convolution kernel W is typically well-localized on 
the sky, which means that the convolution can be com- 
puted directly with a reasonable effort. Alternatively, the 
convolution can be performed by expanding the sky-map 
in spherical harmonics, multiplying the expansion coeffi- 
cients ai m by We, and transforming back. If necessary, the 
high-pass filtering can be further accelerated by making a 
flat-sky approximation locally and performing the convo- 



lution using two-dimensional fast Fourier transforms. All 
these cases involve performing sums rather than integrals, 
since real-world maps are discretized into pixels of finite 
size. This discreteness should not substantially degrade 
the foreground removal as long as the map is properly 
oversampled: Figure 2 shows that the filter has no power 
below the beam size, where C" mse blows up, and therefore 
will not vary appreciable from one pixel to the next. 

We found that this simple method can enhance the de- 
tectability of point sources by a substantial factor, which 
mainly depends on the available signal-to-noise in the un- 
filtered map. The detectability can of course be further 
improved by taking an appropriate linear combination of 
maps at different frequencies (as in Tegmark 1998) before 
the filtering, to remove e.g. the CMB at that stage at the 
cost of a slight noise increase. 

Our filtering procedure has an additional advantage. 
Since the outlier removal scheme is more likely to throw 
away pixels in CMB hot spots than in cold spots, it in- 
evitably introduces a slight bias, correlating false positives 
and false negatives with the CMB. With a 5 — a threshold, 
a 3 — a point source will get removed if it resides in a 2 — a 
CMB hot spot, whereas it takes a 7 — a point source to 
be detected in a corresponding CMB cold spot. By elim- 
inating most of the CMB before the point source removal 
step, our filtering scheme also eliminates most of this bias. 
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FIG. 2 — Same as Figure 1, but in the Fourier (multipole) do- 
main, is plotted (shaded) together with the total power spec- 
trum Cp 0t (heavy line) and the various components that make it up. 
Everything is for Planck channel 5, i.e., at 100 GHz. 

4.1. Robustness of the results 

How sensitive are our results to the various assump- 
tions that we have made? We found that the Galac- 
tic foreground model made only a minimal difference for 
most Planck channels, where the contribution to C\ ot (and 
hence to the confusion noise a) is likely to be dominated 
by CMB and detector noise for all I. Likewise, we found 
that changing the cosmological model had only a minor 
(< 20%) effect on a. MAP and Planck can of course use 
the model determined by their data, iteratively. Instead, 
almost all the uncertainty comes from the assumed source 
count model (To98), so let us compute this dependence ex- 
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plicitly. If we can remove all sources brighter than a flux 
cut S c , the point source power spectrum in the original 
(unfiltcred) map is (TE96; Tegmark & Villumsen 1997) 




independent of £, where n(S) is the number of point 
sources per steradian with a flux exceeding S and [3 = 
—d\nn'/dhiS is the logarithmic slope of the differen- 
tial source counts n' = —dn/dS. This neglects point 
source clustering, which is generally a good approxima- 
tion (To98). The rms fluctuations a ps are given by the 
familiar expression 

i v 7 

where the sum can be accurately approximated by an in- 
tegral. Above we saw that outlier removal gave S c = 
vo jcB(X) = 2-KQ 2 va /c. Substituting this and equation (JsJ) 
into equation (^), we obtain the useful result 



°vs « \j— p Nl ' 2 ^ ( 10 ) 

where N = ird 2 n(S c ) is the number of sources removed 
per beam area. Since relevant values for (3 are typically in 
the range 1.5-2.5 (see references in Tegmark & Villumsen 
1997), the first term is of order unity. Table 1 showed that 
the best attainable a was typically 3-5 times the rms noise 
<r n . Point sources have only a minor impact on a CMB 



experiment if a ps <C <J n , since their power spectra have 
the same shape as that for detector noise (apart from the 
noise blowup up the beam scale) . Equation ( |l0| ) therefore 
tells us that using the CMB map itself for point source re- 
moval is quite adequate as long as N -C (4 x 5)~ 2 = 0.002. 
Conversely, if there are more sources per beam than this 
rule of thumb indicates, then an external point source tem- 
plate will be needed to reduce the point source contribu- 
tion to a subdominant level. This criterion thus partitions 
CMB experiments into two classes: those for which inter- 
nal cleaning suffices and those which need external point 
source data. 

Since a ps tx A 1 / 2 , a useful measure is the safety mar- 
gin, defined as M = (cr n /a ps ) 2 . This is the factor by which 
the number of point sources can be increased before they 
dominate the noise rms. In Table 1, it is seen to be of or- 
der 10 for channels 7 and 8, which means that the models 
of To98 would need to be off by an order of magnitude 
for point sources to imperil the Planck mission. For those 
channels where M < 1 (a ps > a n ), it will be desirable to 
use a multi-frequency subtraction scheme or external point 
source catalogs to further reduce a ps . Future radio source 
surveys at CMB frequencies will therefore be extremely 
valuable to the CMB community. 
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Table 1 — The rows show how the 10 Planck channels can be used for point source detection before (B) and after (A) filtering. The 
values of rms confusion noise correspond to detector noise (cr n ), Galactic foregrounds (o" Ga i is the combined contribution of dust, free-free and 
synchrotron emission), CMB (o"cmb)i the quadrature sum of all of the above (a) and unremoved point sources (<T ps ). <tcmb is computed from 
the full power spectrum - using o"cmb= 10~ 5 T cm ;, ~ 27^tK as in e.g. To98 gives values about 5 times too low. The gain is the factor by which 
filtering reduces cr. N(> 5cr) is the number of detected sources in 8 sr (having S > S c ) and M the factor by which the To98 source counts 
would have to be increased to 
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